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Abstract 

In this paper, we analyse a model for the growth of three-dimensional 
walled cells. In this model the biomechanical expansion of the cell is 
coupled with the geometry of its wall. We consider that the density of 
building material depends on the curvature of the cell wall, thus yield¬ 
ing possible anisotropic growth. The dynamics of the axisymmetric 
cell wall is described by a system of nonlinear PDE including a nonlin¬ 
ear convection-diffusion equation coupled with a Poisson equation. We 
develop the linear stability analysis of the spherical symmetric config¬ 
uration in expansion. We identify three critical parameters that play a 
role in the possible instability of the radially symmetric shape, namely 
the degree of nonlinearity of the coupling, the effective diffusion of the 
building material, and the Poisson’s ratio of the cell wall. We also 
investigate numerically pattern selection in the nonlinear regime. All 
the results are also obtained for a simpler, but similar, two-dimensional 
model. 

Keywords : cell growth modeling, stability analysis of PDE, numerical 
scheme for PDE resolution 


1 Introduction 

The physical features that account for the acquisition and maintenance of 
cellular shape is a current problem in experimental and theoretical biology 

IS El US]. 

In this work, we investigate a generic biomechanical model for the growth 
of walled cells, such as plant cells, fungal hyphae, or fission yeast {S. pombe). 
Cell wall can be described as a thin shell subject to a high internal pressure, 


1 


called the turgor pressure, that can reach up to 10 atmospheres [TB]. Due 
to this mechanical feature, the cell wall by itself determines the shape of the 
cell. It is commonly admitted that a good representation of an expanding 
cell wall is provided by an inflating balloon [H El- 

Several attempts have been made recently to understand the dynamics 
of growth of walled cells, and in particular rod-like cells and pollen tubes, 
from a physical and geometrical viewpoint la El El SO]. Here, we aim to 
further analyze the coupling between the mechanics of cell wall expansion, 
and the pattern of growth, the latter being determined by the cell geometry. 
For this purpose we study a minimal model accounting for the mechanics of 
wall expansion, and heterogeneous distribution of growth along the wall as 
a function of its geometry. 


We follow the works of Dumais et al [5] , and Drake and Vavylonis [B] for 
the biomechanical part of the model. They derive the same set of equations, 
based on slightly different hypotheses. In the former, the mechanical defor¬ 
mation of the cell wall is described using viscoplasticity theory (see also [2] 
where it is described as a viscous thin shell). The latter develops a model 
where the wall is viewed as an elastic membrane under pressure, and sub¬ 
ject to local remodeling. Apart from the equations, the common features of 
these models is the heterogeneity of the mechanical characteristics (either 
the viscosity in Dumais et al, or the rate of remodeling in Drake and Vavylo¬ 
nis). This heterogeneity (denoted as 'k in the present article) depends upon 
the distribution of some growth factor {e.g. Cdc42 for the fission yeast, wall 
building polymers for fungal hyphae cells, or cell wall loosening enzymes for 
plant cells). In this work, we do not attempt to describe with much de¬ 
tails the set up of heterogeneity. We shall describe this heterogeneity factor 
with a generic quasi-stationary reaction-diffusion equation with a geometric 
source term. 

We restrict to the case of an axisymmetric cell. This choice obviously 
rules out many possible shapes, but it is compatible with the acquisition and 
maintenance of the rod shape, which is already a challenging problem [3]. 
Furthermore, the analysis, and even the setting of the mathematical model 
is dramatically simplified in this context. The cell boundary is described 
by the angular deviation of the normal vector from the axis of symmetry, 
ip{t, s), or equivalently by the radius of the cell in the transversal direction, 
r{t, s), where s is the curvilinear abscissa (see Figure^. The biomechanical 
model is a system of two equations which determine the velocity of the cell 
wall in the Frenet frame v = (vn, Vr), 


dvr 

ds 


/ Ks C0S(V9) \ 
\K 0 r ) 


UrCOs((i3) , , , 

-4'(fT0 - i^(Ts) = 0 , (1.1) 

r 

Vr + —^ (<70 - v(Ts) {(Tg- vae) = 0 . (1.2) 

Ke 
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Equations are derived from the constitutive laws governing the deformation 
of an elastic cell wall under constant turgor pressure P [B]. The stresses ag 
and ao are calculated so as to balance the turgor pressure P, 

_ P _ P 

^ 26k0 ’ ^ 25n0 V 

where Kg and kq denote the principal curvatures, and 6 is the thickness of the 
wall. As usually done, the latter is assumed to be constant. We aggregate 
in a single parameter ik the (inhomogeneous) mechanical properties of the 
wall, and the local remodeling rate. We will also refer to 'k as the cell wall 
extensibility, in reference to Dumais et al [ 8 ] . Our main modelling hypothesis 
is that 'k depends on the local geometry of the cell, e.g. via its curvature. 
More precisely, we assume that 'k = P(^), where P is a certain nonlinear 
function, and /r represents the distribution of growth signal [B] at the cell 
wall. Equation for // writes 

- 'jAsnit, s) + s) = PK{t, s ), (1.3) 



where K = KgKg is the Gaussian curvature, and a, /3 ,7 are positive con¬ 
stants. 

We assume that the dynamics of release occurs faster than growth, 
thereby the reaction-diffusion equation (1.3) is at quasi-stationary equilib¬ 
rium. The choice of the source term j3K is motivated as follows: 1) the 
material is released preferentially in the regions of higher curvature, 2) K 
is an intrinsic invariant of the surface, 3) the whole quantity K is con¬ 
stant. Thus, /i(t, s) can be viewed as the result of redistributing a limited 
amount of growth material according to the local geometry, together with 
lateral diffusion on the surface. 

Behind this particular choice for the source term, we have in mind the 
following more realistic process, involving the distribution of microtubules in 
fission yeast (see for instance the computational model developed in 121). It 
is well-established that the cytoskeleton controls cell polarity and cell shape 
HEI]. In normal conditions, microtubules align preferentially in the axis 
of growth. Thus, they deliver a group of proteins (the -|-TIP complex) to 
the cell tip. This enhances the local organization of the growth machinery 
toward cell tips (see m and references therein). Interestingly enough, it is 
possible to redirect the location of the growth zone by mechanically acting 
on the cytoskeleton [m Ei]. These experiments show evidence of a feed¬ 
back loop between cytoskeleton organization and the cell shape. Our set 
of equations can be viewed as a minimal model accounting for this loop. 
We ignore secondary feedbacks that regulate the size of the polar cap. We 
simply assume that microtubules accumulate and deliver the growth ma¬ 
terial in the region of higher curvature. The growth material is retained 
in the polarisome. Consequently, growth is focused at the cell tip. In our 
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basic model, this retention process is associated with the length scale of the 
diffusion process, namely ( 7 / 0 ;)^/^. 

Previous works studied the maintenance of rod-like cell shape. Gener¬ 
ally, the growth pattern is prescribed as a function of the distance to the 
growing tip of the cell [2], e.g. via a Gaussian distribution of the growth 
material along the cell wall. In p], the variance of the Gaussian distribu¬ 
tion can depend on macroscopic quantities such as the length and the mean 
radius of the cell. Here, on the contrary, we do not assume that the cell 
is rod-shaped initially. Moreover, the distribution of the growth material 
intrinsically depends on the local geometry of the cell wall, e.g. through 
its curvature. We focus on the possible initiation of a rod shape from a 
spherical one, as observed experimentally mi- We address mathematically 
the following morphogenesis question: starting from a small perturbation of 
a growing radially symmetrical shape, can the model evolves towards a rod 
shape? We fully answer this question using refined linear stability analysis 
of the radially symmetrical growing shape. 

Other studies have investigated the dynamics of cell growth in prokary¬ 
otic actinomycetes. On the contrary to eukaryotic fungi, such as fission 
yeast, the cytoskeleton plays certainly a minor role in the establishment of 
the growth pattern: wall building material is likely to be transported to 
the tip by diffusion. Goriely and Tabor propose two models to investigate 
self-similar growth of the tip. In mm, they develop a model based on 
large-deformation elasticity theory. The membrane is described as an ax- 
isymmetric elastic shell far from the reference configuration. Similarly to the 
models discussed above PE], the inhomogenous elastic modulus is given a 
priori, being assumed that the wall gets stiffer far from the tip. In uni, 
the same authors develop a purely geometrical model for tip growth, thus 
neglecting biomechanical effects. The models follows previous studies on the 
morphogenesis of unicellular algae IlsiDrai. In this model, surface evolution 
is determined by purely kinematical considerations, and local deposit of ma¬ 
terial, yielding areal growth. The rate of areal growth is a function of the 
local geometry {e.g. the Gaussian curvature). The authors obtain various 
shapes of self-similar growth in the 2D and in the 3D case. 

We can summarize our results as follows. Firstly, we noticed that without 
taking into account the mechanical effects, the dynamics of wall expansion 
in the normal direction proportionaly to [i (or F{p)) instantaneously creates 
a very singular pattern. In fact, the system is likely to be ill-posed. Unsur¬ 
prisingly, taking into account mechanics has a stabilizing effect. We are able 
to measure quantitatively this effect in the linear regime. The range of pa¬ 
rameters for which the system undergoes symmetry breaking is surprisingly 
narrow. For instance, it is required that F is strongly nonlinear. In addi¬ 
tion, 7 , the diffusion coefficient of fj,, has to be relatively small. Then, the 
morphogenetic instability arises as a competition between the mechanical 
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effects, and the redistribution of material according to the geometry. 

Our analytical work is complemented with numerical simulations far 
from the radially symmetric shape. Interestingly enough, it is unstable, the 
system selects various anisotropic shapes, including the rod-like shape. 

We believe this work paves the way for further (nonlinear) investigations 
of the coupling between growth and form, and provides some quantitative 
basis for the analysis of shape selection in walled cells. An interesting per¬ 
spective would be to adapt this framework to other types of growing walled 
cells, such as bacteria E. coli [3], for which another kind of feedback between 
cell wall growth and curvature has recently been established [23]. 


The paper is organized as follows. Section deals with the complete 
presentation of the model. After introducing some geometrical notations in 
Subsection 12.11 kinematics are described in Subsection 12.21 and the mechan¬ 


ical aspects are given in Subsection 2.3 The last Subsection 2.5 introduces 


an analogous model in the two-dimensional case. It governs the growth 
of a closed curve, including the main ingredients of the three-dimensional 
axisymmetric case. Section presents the main theoretical results of this 
paper. We deal with the 2D case in Subsection 3.1 The linear stability 


analysis of the three-dimensional case is expressed in Subsection 3.2, In ad¬ 


dition, under the (unrealistic) hypothesis of a fixed length, we characterize 
in Subsection |3.3| the stability of the spherical configurations. The proofs 
are given in Section Section presents numerical results. We confirm 
the results of the linear stability analysis, and we illustrate the nonlinear 
dynamics of the system. 


2 Description of the biomechanical model 

In this part, we introduce step by step the PDE system which governs the 
evolution of the cell wall. We focus on a model which combines two main 
conceptual parts. 

1. The cell wall expansion derives from the constitutive laws which govern 
the deformation of an elastic boundary under pressure, and subject to 
remodeling. We use a standard mechanical model proposed in mw- 
The mechanical parameters are the wall extensibility 'k, the Poisson’s 
ratio and the turgor pressure P. 

2. The extensibility depends on the local concentration n of some growth 
material which is distributed along the cell wall. The distribution of 
material depends intrinsically on the geometry of the cell. We opt for 
a basic coupling via the Gaussian curvature as in uni, but the math¬ 
ematical results obtained here can be generalized to other couplings. 
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2.1 Geometry of the cell wall 

We restrict to axisymmetric cells. The cell wall is a surface of revolu¬ 
tion. The generatrix curve in the x-z plane is parametrized by Ct := 
{{r{t, s), z{t, s))} (see Figure where s is the curvilinear abscissa. We de¬ 
note by L(t) the length of the generatrix curve at time t > 0. As described 
in Figure ip is the angle between the normal of the curve n and the ^;-axis, 
while T is the tangent vector. We note (61,62,63) the canonical basis of M^. 
We also denote by Ks and Kg the two principal curvatures of the surface. 
Since St is a surface of revolution, they have the following expression (see 

m p- 161 ), 


_ - (5^r)2 

r{t,s) 


Kg{t,s) 


-dsv 

(5^r)2 


( 2 . 1 ) 



Figure 1: Generatrix curve Ct of the surface of revolution in the a; 2 :-plane 


2.2 Cell wall expansion: kinematics 

For the sake of completeness, we derive the equation governing the evolution 
of a curve which moves according to a vector field, v(t, •) : G —)■ M^. 
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Figure 2: Infinitesimal element, called ds, pushed by a vector field v. 


Let s £ [0,L{t)]. According to Figure]^ we express the time variation 
I (|§) as follows, 

d^r{t,s) dr{t,s) fdv{t,s) \ dv{t,s) 

■ 

Moreover, by following the variations around a material point, we get 

d /dr{t, s{t))\ d‘^r{t,s{t)) ds{t) d‘^r{t, s{t)) 

dt \ ds ) dtds dt ds"^ 


On the other hand, we notice that the curvilinear abscissa is reparamatrized 
by the elongation of the curve. 


ds{t) 9v(t, s') 


dt 


ds 


T{t, s') ds', 


Finally, the equation governing the evolution of the generatrix curve reads, 
for s £ [0,L(f)], as follows 


d fdr{t,s)'\ Off f^d^r{t,s') /^ j A ^ , ^^r{t,s) 

m[-d^) = TA~yi + 

( 2 . 2 ) 


We refer to Appendix [A| for more details concerning the derivation of (2.2). 
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2.3 Cell wall expansion: mechanics 

The cell wall dynamics derives from the constitutive laws governing the de¬ 
formation of an elastic boundary under pressure. We refer to the biophysical 
literature for the justification and more in-depth derivation of the equations 

jZlEllElES]. 

Following [S], we consider the following constitutive equations relating 
the strain rates and the stresses, 

I e. a,) , 

\ 60 = {ae - as) , 


where ag (resp. ae) is the meridional stress (resp. the circumferential stress). 
Parameters are: v, which is analogous to the Poisson ratio in linear elasticity, 
and dt, the cell wall extensibility. We refer to [ 6 ] which derive an equivalent 
model to the one contained in [5] in a slightly different context (elastic cell 
wall remodeled under turgor pressure). In [Fj, 'h is linked to the remodeling 
rate of the cell wall. 

The assumption of axial symmetry enables to relate the stresses ag and 
a 0 to the principal curvatures Kg and kq 


ag 


<70 


P 

2 k 0 ’ 
2k0 V 



(2.4) 


Kinematic relations for axisymmetric shells make the relation between the 
strain rates and the velocity. As usual, we decompose the velocity field v as 
V = u„n -|- VrT. We have [B] 


ee 


— T 


= Vnl^e -h 


dVr 
ds ’ 

Vr COs((/?) 
r 


(2.5) 


Thus, by substituting equation (2.3) into (2.5) and by using the relations 


(2.41, we deduce that the functions 
equations on M x [0,L(t)], 


and Vr solve the following system of 


dvr 

ds 


Vr cos(ia) 

i^eVn H- ^{cre - I'CTg) = 0, 

r 

(Kg cos((p)\ Kg ^ , s ^ , N „ 

- Vr H-^ (<70 - uag) -^{ag- vae) = 0 . 

\hi 0 r J Ke 


( 2 . 6 a) 

( 2 . 6 b) 


2.4 Cell wall expansion: growth pattern 

We assume that the function \k, which can be read as the cell wall exten¬ 
sibility, or the remodeling rate, depends on some growth factor released in 
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the vicinity of the cell wall. The surface concentration of this growth factor 
is denoted by s). We assume a general relationship between 'I' and //, 


where F is an increasing function. We suppose that the growth factor ^ 
is diffusing laterally, and is released locally, proportionally to the Gaussian 
curvature. Moreover, we assume that the dynamics of release-diffusion are 
much faster than the time scale of remodeling and growth. Thus, the density 
/i satisfies the following equation 

- 7 Asn{t, s) -h a n{t, s) = (3Ks{t, s)K 0 {t, s ), (2.7) 


where 7 ,a and 13 are positive constants and As is the Laplacian-Beltrami 
operator on the cell wall S. In the case of an axisymmetric surface, the 
Laplacian-Beltrami operator As of ^ is expressed as follows lap. 79]. 


- As^i{t,s) = - 


r{t,s) 




S) 

ds 


( 2 . 8 ) 


At this stage, we can compute the rate of expansion of a spherical shape. 


Proposition 2.1 (Spherical solution) Let L'^ be the solution of the fol¬ 
lowing ODE 


5 ^ = ^p(i-OF 

dt 2 ^ ^ 


/^TT^ 


L%tf 


with a given initial length L'^(O). Then, defined by 


= 


TT 


/Svr^ 


ct+\2 ’ 


aL^{t) 


is a partieular solution of the system ( 2 . 2 |-( 2 . 6 )-( 2 . 8 ). 


(2.9) 


( 2 . 10 ) 


In Section]^ we address the issue of stability of the solution (2.10). 


2.5 A simplified 2D model 

In this subsection, we derive the dynamics which governs a thin elastic string 
under pressure. More precisely, we assume that the cell wall is a closed curve 
which is symmetric with respect to the z-axis. 

We first express the Hooke’s law, e = ikcr, where e is the strain rate 
and a is the stress acting on the curve. The latter is determined by the 
Laplace law, a = P/k, where k is the curvature. Moreover, by geometrical 
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consideration, the strain rate is also equal to e = where Vn is the normal 
component of the velocity vector field. All in all, the normal velocity of the 
curve is given by, 

= ( 2 . 11 ) 

We use the same notations as for the 3D axysimmetric case. In particular, 
we have k = —dg^p. By analogy with the 3D case, the evolution of the curve 
obeys the following equation 


d /'dr{t,s)\ 
(9t V ds ) 


ds 



dgifit, s' 


■ds' 


dr{t, s) 


ds 


+ 


^ (p 


ds\ {ds(p{t,s)y 


sm{ip{t,s)) 


( 2 . 12 ) 


The coefficient 'h is supposed to depend on some growth material which 
is released proportionally to the curvature, and diffuses quickly along the 
curve. 


5V(i,s) adip{t,s) 

7 — ^79 — + 'S) = 


(2.13) 


Qg2 ' t' 

The next proposition states that the circular shape is preserved by the 
system (2.12)-(2.13|. 


Proposition 2.2 (Circular solution) Let L'' be the solution of the fol¬ 
lowing ODE 


dL%t) 

dt 



L%tf 


(2.14) 


with a given initial length T'^(O). Then, {r^,pL) defined by 


, L''(t) . f STT \ 

r(M) =—™(^) 




TT 

fj-K 


aL^{t) ’ 

is a solution to the system of equations ( |2T2| )-( |2T3t . 


3 Linear stability analysis 

In this section, we derive the precise conditions under which the cell wall 
expansion is unstable around the radially symmetric solution. The result 
is given in Subsection 3.1 for the 2D case and in Subsection 3.2 for the 3D 
case. Note that the 2D model cannot be viewed as a particular case of the 
3D system, but there are strong similarities between the two settings. 

To analyze the stability of the model, we proceed in two steps. The 
first step consists in transforming the original system into a more tractable 
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one by using an appropriate change of variables and by making the physical 
quantities dimensionless. The second step is to solve the linearized equation 
around the circular (2D case) or the spherical (3D case) solution. 

The following definition will be useful. 


Definition 3.1 (Degree of nonlinearity) For f : M+ —)■ M a differen¬ 
tiable function, and fa E M+, we define the degree of nonlinearity deg(F; ^) 
of f at point p, as following 


deg(F; fi) 


f{F)F 

fiF) 


dlogfil^) 

dlogifi) 


It coincides with the notion of elasticity in economics. For a homogeneous 
function f{fj.) = pP, we have deg(F; /i) = p. 

In the sequel, the curve is parametrized by the angle (p. The following 
relation enables to rewrite the equations for the wall expansion (2.2) and 

dr{t, s) 


(| 2 T 2 l), 


ds 


= cos{(p{t, s)) 


(3.1) 


3.1 Stability results for the 2D model 

Firstly, we adimensionalize the equations. Let p and fi be defined as 

p{t,x) = p{t,xL{t)) and /2(t, x) = xL(t)) =, t>0. 

We introduce ^(t, x) = F (^^f^iF{t)il{t, x)^. We also introduce the reduced 
diffusion coefficient on the membrane. 


a{t) 


aL{ty 


Then, {(p, fi) is solution to the following system of equations 




^(t,u) 


L(t) 

- pm 


0 dx<p{t,u) 

dx^{t,x) 2§(t,x) 2 


{dxp{t,x)) {dx<p{t,x)y 


d^p{t,x) 


(3.2) 


2 dlfL{t, x) + jl{t, x) = Pdxp{t, x ), (3.3) 


cr(t)„o_,. . ..... 1 

TT^ 

with the boundary conditions dxfi{t, 0) = dxft{t, 1) = 0. The following result 
quantifies the local stability of the system (3.2)-(|3^ around the circular 
solution. 


X E (0,1). 
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Theorem 3.1 (Linear stability analysis) Let {u,v) be the solution of 
the linearized version of system (3.2 )-( [3^ around the eireular solution. 
Let {ak)k>i denote the Fourier eoeffieients of dxU, defined as dxu{t,x) = 
akjt) cos(fc7rx). The eoeffieients are solutions to the following set of lin- 

k>l 

ear ODE 

a'kit) = M{t,k)ak{t) , 

where M is a rational fraetion of k, the sign of whieh is determined by the 
sign of the following polynomial for any integer value k > 1, 


N{t,k)= I - deg{F-, n^{t)) + cr^{t) - 2 + deg{F-, n^{t)) 


k^- 


2a^{t) 


k^ 


The proof of this result is given in Subsection 4.1 We define the functions 
Gi and G 2 as, 

Gi(cr, d) = 1 + (d — 2) — 

a 

02 ( 0 ', d) = 1 + 2 (2 — 3d) —h (d — 2) —^ . 


a 




For a given t > 0, the zeros of N{t, ■) are expressed as 


1 

zb — 

2 


^Gi (cr':(t),deg(F;/ryt))) ± ^G 2 {a^{t),deg{F-, n'^{t))). 


The consequences of this result are discussed in Section 3^ and illustrated 
in Section [H 


3.2 Stability results for the 3D model 

In this Subsection, the previous results are extended to the 3D case. We 
proceed with the same nondimensionalization. Let (p and jl be defined as 

<p{t,x) = ip{t,xL{t)) and fi{t,x) = p{t,xL{t)) = ’ t>0, x G (0,1) . 

We introduce x) = F We also introduce the reduced 

diffusion coefficient on the membrane. 


a{t) 


aL{ty ’ 


and the auxiliary functions Icos and Icosin, defined as 


rx 

Icos(t, x)= / cos{if{t,x'))dx' and Icosin(t,x) 

Jo 


Icos(t, x) 
sin((^(t, x)) 


(3.4) 
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Then, {ip{t,x),jl{t,x)) is solution to the system of equations 


dtif{t,x) = {d{t,x) + H{t,x)dj:<p{t,x)) dx‘f{t,x) 

+ cot((^(t, x)) d{t, x) — dxH{t, x) , (3.5) 


fa2~ , cos((^) \ 1 ,sin((^) 


vr^ 


dxi^ 


Icos 


dxp' ) “1“ A — 9 dx^ 


7T^ 


Icos 


(3.6) 


where the function d and H are given by 


1 


1^(1 — 2i/)Icosin + 2 (i/ — 1) (Icosin)^ dx^f + (Icosin)^ (dx(p)^'j , 

(3.7) 

"" " ^(2 — z^) (Icosin)^ — (Icosin)^. (3.8) 


Ai = -FL^i(l- 


Ao = -PLT (2 - 


The details of the change of variables are given in Appendix The 


following result quantifies the local stability of the system (3.5 )-(3.61 around 
the spherical solution. 

Theorem 3.2 (Linear stability analysis) Let (u,v) be the solution of 


the linearized version of system (3.51 -(13.6^ around the spherieal solution. 


Let {ak)k>i denote the Fourier eoeffieients of dxU, defined as dxu{t,x) = 
akjt) cos(fc7rx). The subset of eoeffieients A{t) = {ak{t))k >2 solu- 


k>l 


tions to a triangular system of ODE, denoted by 

Aft) = M{t)A{t) , 


(3.9) 


where the signs of the diagonal eoeffieients of M are determined by the sign 
of the following polynomial for any integer value k >2, 

N{t,k) = nof) + ni{t)k + n 2 {t)k‘^ + n 3 {t)k^ + n 4 {t)k^ , Vt > 0 , (3.10) 

with 


no{t) = (1 - n) [1 - 2deg(T; M^i))] , 
ni{t) = -(1 - n ) [a^{t) + deg(F; fi'^it))] + 1 , 
n 2 {t) = -va^{t) - 1 + (1 - n) deg(F; pA{f )), 
n^{t) = 2cr''(t) , 

= —cr^{t). 
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The proof of this result is given in Subsection 4.2, We define the functions 
Gi and G 2 as 


Gi(a, d)=2,-2v + 2{{l-v)d-l)-, 

a 


1 


1 


G 2 (cr, d) = {1 - v) - 2(1 - v){{2> + v)d - 1)—h ((1 - v)d - 1) ^ . 


a 




For a given t > 0, the zeros of N{t^ ■) are expressed as 


1 1 

— i — 

2 2 




3.3 Instability in the case of quasi-stationnary expansion 

In this Subsection, we discuss several necessary conditions for instability 
based on the previous analysis. 

We make two assumptions to simplify the discussion. Firstly, we do 
not restrict k to taking integer values. We assume that it can be any real 
value greater than unity. Secondly, we make the modelling assumption that 
the dynamics of (in)stability is faster than wall expansion. Therefore, we 
consider that L(t) varies slowly, so slowly that we further assume that it 
takes a constant value L. 


The 2D case. We begin with the two-dimensional case. There are four 
possible roots, which can be real or not by pairs. In addition, we require 
that the largest root is greater than one. 

We find that the latter condition implies Gi + ^/G 2 > d, which is satisfied 
if Gi > 4 or if G 2 > (4 — Gi)^. The latter cannot be satisfied since G 2 — 
(4 — Gi)^ = —8(1 -|- l/u) < 0, therefore we obtain the necessary condition 
Gi > 4, which reads 

d > 2 -|- 3(T . 

In particular, it is mandatory that the degree of nonlinearity is larger than 2. 
This could have been guessed from the expression of the expansion velocity 
. In fact, if 'k = then the velocity of expansion is constant, and 
the geometric coupling disappears. 

On the other hand, a simple asymptotic analysis of the roots in the 
two cases i) d —>■ -|-oo, and ii) cr —)■ 0 together with the conditon (d > 2), 
shows that high degree of nonlinearity and small lateral diffusion on the 
cell wall guarantee instability. In fact, by introducing the notation ki = 

yGi + y/G 2 and k 2 = yGi-V^, we notice that i) as d —>■ -|-oo, we have 
respectively ki —)■ -|-oo and A ;2 —^ 1, and ii) as cr —)■ 0, we have respectively 

ki —)■ - 1-00 and k 2 ^ > 1. 
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The 3D case. We perform the same analysis in the three-dimensional 
case. Notations are the same as in the 2D case. We find that the condition 
ki > 2 implies Gi -|- 2^/G2 > 9, which is satisfied if Gi > 9 or if 46*2 > 
(9 — Gi)^. The latter cannot be satisfied since 4^2 — (9 — Gi)^ = —16(1 -|- 
i>){2 + 1 /a) < 0. Therefore, we obtain the necessary condition Gi > 9, 


which reads also 


d> 


1 -|- (3 -|- l')(T 
1 - u 


(3.11) 


The same asymptotic analysis of the roots holds true, provided {l — v)d— 
1 > 0, which is compatible with the above necessary condition. In particular, 
i) as d —)• -|-oo, we have respectively ki —)■ -|-oo and k 2 —)■ 2, and ii) as a —)■ 0, 
we have respectively ki —)■ -l-oo and k 2 ^ \ + > 2. 


4 Proofs of the linear stability results 


This Section is devoted to the proofs of Theorem 3.1 and Theorem 3.2 


respectively. Although the three-dimensional case requires extra compu¬ 
tational tricks, the proofs are sketchily the same. We drop the ~ in the 
adimensional system for the sake of notation. 

In the following, we denote respectively by /r"', and small per¬ 
turbations of the radially symmetric shape associated whith the particular 
solution (p{t,x) = TTx, iJ,{t,x) = 1, and 

We use the Fourier basis associated with periodic functions of period 
two, up to extending the functions for x G (0,2) using axisymmetry, e.g. 

^p{x) = 277 — (p(2 — x) and fi(x) = /r(2 — x), Vx G (0,2). (4.1) 

4.1 Proof of Theorem 13.11 : 2D case 


We aim at linearizing the system (3.2|-(3.3| around the circular shape solu¬ 
tion (ifjfJ.) = (vrx, 1). The linearized equation of (3.2) is 

^"'dx' + 


TT^ 

\ vr JO 

xdxP 


I m) 


X TT 


mJ 


(4.2) 


Since L'^ solves equation (|2.14|), equation (|4.3|) becomes 

^dxdf 


dtp^ = -PL^{t) l/^'dx' + ^dx 


(4.3) 
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From the definition of in the nondimensional variables (see Section 


3.1), we have 




L-{t) 


= deg{F■,^x^{t))F{^^%t))[^r{yx)- 


(4.4) 


By differentiating (4.3) with respect to x, we get 


dt x)) = - PL''(t)deg(F; ^f{t))F{if{t)) /i"’(t, x) + x) 


TT^ 


+ PL%t)F{^,%t)) f-a,<^-(t,x) + ‘^dlip^{t,x) 

\'K 'K'^ 

+ Pdeg{F■,^i%t))F(^x%t))L^{t). 


+ vr 


(L'{t)V 

\L{t)) 


On the other hand, equation 
ayt) 


TT^ 


is already linear, so we write 
1 

— I 

TT 


dliF’{t,x) + fF’{t,x) = -dx^'^{t,x) 


(4.5) 


(4.6) 


We expand the 2-periodic functions dx^'^ and [jF in Fourier series. Using 
the symmetric relations (4.1), we can write 

dx^p^{t, x) = ^ ak{t) cos{kT:x ), 
fc>i 

= Y, Ck{t) cos(/i:7rx), 


k>l 


where 


ak{t) = 2 / cos(A;7rx)9a;¥^’^(i, 2:)dx, 

Jo 

Ck{t) = 2 / cos(A:7rx)/i“’(t, x)dx. 

Jo 

By replacing the functions dx^'^ and fjF by their Fourier series in equation 

1 


(|4.6|) we get 


Ck{t) = 


TT (1 + (T’^{t)kP‘) 
Similarly, for equation (|4.5|), we obtain. 


dak 1 


ak{t), V/c > 1. 


(4.7) 


(-7rdeg(F; //<=) (l - k^) + (l - 2k^) a^) . 

Therefore, we deduce that for all A: > 1, the following differential equation 
holds, 

dak PL'^F 


dt vr (1 + a^k"^) 
which concludes the proof. 


((1 - deg(F; + (rx'' - 2 + deg(F; k^ - 2a‘'k^'j ak , 


L^{t) 
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4.2 Proof of Theorem 13.21 : 3D case 


We aim at linearizing the system (3.5|-(3.6) around the spherical shape 
solution (y?,/i) = (ttx,!). 

Step 7^1 : Linearized equations. We start by expressing the expansion 


of functions d (3.71 and H (3.8) around the spherical solution. The zeroth- 


order term reads 




The first-order term is 


= (1 - z.) ^ 




( 21 .- 1 ) 


/ cos (y?'") So^V^^'dx- 

sm((^'^j Jo vr 


Similarly, the expansion of d writes 


k\ = 0 , 


A)" =z/-L'=^" 


TT [ sin {ip^) 


[ cos 

Jo vr 


Then, after differentiating the linearized version of equation (3.5) with re¬ 
spect to X, we get 

dt (dxpn = ((x)" " ^ 


TT 


sin((^' 


fi \2 


A)^+cot(^'=)9,Ar-W. 


The equivalent of (4.4) in the 3D case reads 


x) = deg (F; if{t)) F{if{t)) x) - 2 


L-{t) 


(4.8) 


(4.9) 


We obtain, 

dt {dxif-^) =PL^F{pk) 
1 


sin((^'^)” Jo 


)3. 

2 


T — [1 — r' — — r). . 

27r V sin2((^‘=) 


cos {p^) dxP^ dx 
cos{ip^ 


dxjF 


27r2sin(<^'=)^*^'^ ^ 27r3^^^'^) 
- PL=deg(F; p^)F Q (1 - z.)^ (1 - z.) 

/ rz\ 1 

+ vr(^-j +2(l-i^)PF(2deg(F;//‘=)-l)L“'. (4.10) 
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Then, we multiply both side of equation (4.101 by sin((/?'^)'^ and we differ¬ 


entiate with respect to x. Eventually, we divide the resulting equation on 
both side by sin((/7'^)^, so as to obtain 


dt (Stt cos{ip‘^)dx<f'^ + 


= PL^F\^{1 - v) cos{ip^) dx<p^ - ^{2 + 1 ^) sm{ip'^)dl(p^ 

- PL'=deg(F; - v) |^^7rcos((/p'=)/r“' + ^ sin(<^'=)5,;U- 

+ -{1 - i/)PF(2deg(F;/r^) - 1) 7rcos((^^)L’^ -h Svr^ cos(¥?'^) f — j . 

(4.11) 


By multiplying the linearized version of equation (3.61 by sin((^'^) and by 
differentiating both side with respect to x, we get 


1 


2 cos{(p^)dxip^ -I- — sm{ip‘^)dlip'^ 

TT 


' C^a2..w 


= vr cos(-I-(cr'^-1-1) sm{ip^)dxfJ-'^—2— cos{ip^)dxfJ-^ — ^ sm{ip^)dxfJ-^ 


TT 


TT^ 


(4.12) 


Step 7^2: Fourier series expansion. We expand the 2-periodic functions 


dx<f^ and fjF in Fourier series. Using the symmetric relations (4.1), the 


Fourier series read as in the 2D case. We will need the following formulas: 

cos{ip'^)dx(p^ = ai -|- (ao -|- 02 ) cos(7rx) 

-f (ofe+i -h ak-i) cos(/c7rx) . (4.13) 

k>2 

sm{(p^)dx‘p^ = —vrai — 27ra2 cos(7rx) 

- vr ^ ((/c -h 1) Ofc+i - (/c - 1) afe_i)cos(/c7rx) . (4.14) 

k>2 


cos{(p‘^)dx(p^ = —TT^ai — 47r^a2 cos(7rx) 

— TT^ ({k -|- 1)^ gfc-i-i -|- (fc — 1)^ Qfc-i) cos(fc7rx) (4.15) 


k>2 
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= vr^ai + 87r^a2 cos(7rx) 

+ vr^ ^ (^{k + 1)^ Ofc+i - {k - 1)^ cos(/c7rx) . (4.16) 

k>2 

In addition, the Fourier series of the functions cos{(p^), sm(ip^)dxfJ-^, 
cos{(p^)dxH'^ and sin((/?'^)9^/i"' are involved in the linearized equations ( |4.11[ ) 
and ( 4.12[ ). They are obtained by replacing the function by for 

k = 1, - ■ ■ ,4 and the coefficients ai by Ci for i > 0, in formulas (4.13)-(4.16). 
We obtain that the Fourier coefficients for k >2 are linked via the following 
relation 

^ {Ml A) = PL^^F {M 2 A + deg(F; p^)MsC) , (4.17) 

where A is the infinite column vector A = {ak)k> 2 : C = {ck)k> 2 : and the 
Mi, i = 1,2,3, are upper triangular infinite matrices 


/pi(2) 0 P2{2) 0 


Ml = 




0 pi{k) 0 P2{k) 0 


V 


(qi{2) 0 92(2) 0 


M 2 = 


0 qi{k) 0 q 2 {k) 0 

(r{2) 0 -r(2) 0 


M3 = 


V 


0 r{k) 0 —r{k) 0 

The polynomials Pi,qi, i = 1,2 and r are expressed for A: > 2 as 


/ 


Pi{k) = TT {2 + k) , P2{k) = TT {2 - k) , 

qi{k) = —k^ — k‘^ + {3 — u)k + 2(1 — u ), q 2 {k) = k^ — — {3 — u) k + 2(1 — v ), 

r{k) = t: {1 — iy)k [k^ — 4^ . 


Since Mi does not depend on time, the set of equations (4.17) becomes 


^A = PL^F 
dt 


(MfiM2A + deg(F;/i‘=)Mf^M3C) . 


(4.18) 


Remark 4.1 (Competition between mechanics and expansion) We no¬ 
tice already some interesting properties. First, the roots of polynomial qi are 
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—2, i ± i-v/5 — 4i/. Thereby, for k >2, qi{k) < 0 sinee v G (0,1). Conse¬ 
quently, all the diagonal eoeffieients of M 2 , whieh are equal to are 

negative. It underlines the faet that the meehanies stabilizes the spherieal 
solution. On the other hand, M 3 represents the eoupling with heterogeneous 
growth. Sinee the roots of polynomial r are —2, 0 and 2, then > 0 

for k >2, heterogeneous growth is elearly a destabilizing phenomena in this 
setting. The present analysis enables to quantify this eompetition. 


By using the linearized equation (4.12), we obtain 


where 


C = N{^N2A, 

(si{2) 0 S 2 { 2 ) 0 


(4.19) 


Ni = 


0 si{k) 0 S 2 {k) 0 


V 


and 


/ti(2) 0 t2{2) 0 


N 2 = 


0 ti{k) 0 t2{k) 0 


V 


/ 


The coefficients Si{k) and ti{k), i = 1,2 are expressed for k>2 as 

si{k) = -K [k — a^k^ + a^k^^ , S 2 {k) = —tt + a'^k^ + a^k^^ , 

ti{k) = k + l, t2{k) = -{k -1) . 


By replacing C in (4.18) by its expression in (4.19), we get eventually 

d 


dt 


-A = M^. 


(4.20) 


where 


M = 


(Mf IM 2 + deg(F; /r")Mf ^MsiVf ^iVs) • (4.21) 


All the matrices appearing in the above expression are in fact upper trian¬ 
gular. The diagonal coefficients of M are denoted by A^. For k > 2, we 
have 


Afc — 




qi{k) + deg{F-, fj.^) 


r{k)ti[k) 

si{k) 


= PL^F 


N{t,k) 
D{t,k) ’ 


pi{k) 

where the numerator is given by 

N{t, k) = no(f) -|- ni{t)k -|- n 2 {f)k^ -|- nz{t)k'^ -|- ni{t)k'^ 
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with 


no{t) = {l-u)[l- 2deg(F; , 

ni{t) = -(1 - u) [a^{t) + deg(F; n‘'{t))] + 1, 
n 2 {t) = - I + {I - v) deg(F; , 

n^{t) = 2cr''(t) , 
rn{t) = . 

and the denominator is given by 

Z)(t, k) = TT {a‘^{t)k^ — a^{t)k + 1^ > 0 . 


5 Numerical results 


The numerical scheme is briefly introduced, and numerical results are dis¬ 
cussed in the 3D case. In particular, we investigate the linear stability of the 
spherical shape, and we explore the behaviour of the system in the nonlinear 
regime. 


We apply standard finite difference methods for transport-diffusion equa¬ 
tions. We opt for a piecewise linear approximation of the function (p on the 
regular grid 0 = xq < xi < • • • < Xm+i = 1 with space step — x* = Ax. 
The functions dxP> and p are approximated by piecewise constant functions 
on the shifted grid Ax/2 = yo < Ui < • • • < Vm = ^ — Ax/2. The trans¬ 
port parts, respectively in (3.2), and (3.5) are approximated by a first-order 
upwind scheme. As for the diffusion terms, we use a semi-implicit approxi¬ 
mation. We have paid much attention to the discretization at the boundary 
in order to remove the possible singularities. For instance, the approxima¬ 
tion of the integral Jq cos (p{t, x') dx' (3.4), defined on the shifted grid, is 
computed using the piecewise linear approximation of p. 

We checked carefully that our numerical results are in agreement with 
the linear stability analysis. 


The initial condition (^(0, x) is a perturbation of the spherical shape 

(p{x) = vrx having ten Fourier modes with random coefficients. 

The coupling function F is homogeneous, F{fi) = so that deg(F; y) = 

d. We shall discuss the influence of d, the degree of nonlinearity, and a the 

membrane diffusion coefficient. An important point is that a depends on t 

2 

due to the increase in length (recall that a{t) = 


Quasi-stationary expansion. Firstly, we suppose that cell expansion is 
slow enough to consider that L(t) is approximately constant. 

Unsurprisingly, the instability region in the {d, a) space lies below a curve 
which is increasing with respect to d (Figure j^a)). Therefore, the instability 
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Figure 3: (a) Stability region in the space {d, a) for u = 0.5. The red curve is the 
set of points such as the difference between the two greatest roots of polynomial 
(3.101 are equal to 1. For the sake of comparison, the necessary condition derived 
in (3.111 is plotted in dashed line, (b) Influence of the Poisson’s ratio ly. 


is favoured as d increases or cr decreases. As a by-product, instability is 
enhanced as the cell is expanding, since a scales as the inverse of L{t)‘^. On 
the other hand, for a given {d,a), the instability region gets larger when u 
increases (Figure [3];b)). 

We present in Figure]^ a numerical simulation for a set of parameters for 
which the spherical shape is stable. We observe that the initial perturbation 
is quickly relaxed towards the constant state dx^ = vr. We present in Fig¬ 
ures 1 ^ two numerical simulations for the same set of unstable parameters, 
but with different initial conditions. This is to illustrate various possible 
nonlinear behaviors. We observed at least two possible profiles, both in the 
two-dimensional case and the three-dimensional case: i) dx^ is above the 
mean for intermediate x (Figure j^a)), or ii) dx'^ is above the mean for ex¬ 
tremal X (Figure j^b)). To explain the selection of two possible modes, we 
opt for a simplified description in 3D. We consider an inflating ellipsoid, the 
axes of which are denoted by a{t),b{t),c{t). By axisymmetry, we assume 
w.l.o.g. that a{t) = b{t). We assume the following rule for the evolution 
of the curve: the rate of expansion of the axes is equal to the Gaussian 
curvature at the corresponding tips. Therefore we write 

da{t) 1 

~ir= 

dcjt) ^ ^ ^ cjtf ^ c{tf 
dt a{tYb{tY fl(0^ 

We notice that the quantity — is constant for t > 0. We denote 
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Figure 4: Evolution of function dx<f for a stable set of parameters cr = 10 ^ and 
d = 4. Initial shape 9a;(/?(0,a;) is plotted in dashed line. 


5 = a(0) ^ — c(0) The first equation in (5.1) becomes 


^-3 


Mt) ^ _ 

dt \a{t)^ ) 


The asymptotic behaviour depends on the sign of 5. If d > 0, i.e. c is the 
major axis of the ellipsoid, then a{t) converges asymptotically towards the 
globally stable equilibrium Ooo = whereas c blows-up in finite time. 

On the contrary, if d < 0, i.e. c is the minor axis of the ellipsoid then a{t) 
grows linearly with time. On the other hand, the second equation in ( |5.1[ ) 
becomes 


dc{t) 

dt 


c(t)2 




therefore c(t) converges asymptotically towards the globally stable equilib¬ 
rium Coo = (—As a conclusion, depending on the relative length 
of the axes initially, the ellipsoid will evolves towards a cigar shape, or a 
flat expanding shape. Of course, all aspects of the mechanics have been 
neglected in this toy model. 


Inflating system. Then, we investigate the evolution of the system in the 
case where L(t) is free to increase (when the rate of expansion is comparable 
to the relaxation time of the system). We expect that the system becomes 
more likely unstable as L{t) increases (see Figure]^. We also expect more 
sophisticated patterns. In Figure we use the same set of parameters and 
initial condition as in|^ We observe the same dynamics for earlier time, but 
the spherical shape becomes progressively unstable, as expected. 

Finally, we present in Figure two possible behaviours in the nonlinear, 
inflating regime, as in FigureWe also show the reconstruction of the three- 
dimensional axisymmetric cell as a result of the numerical simulations. We 
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Figure 5: Evolution of dx‘f for an unstable set of parameters a = 5.10“^ and 
d = 4, and two different initial conditions. Initial shapes are plotted in 

dashed line. 




Figure 6: Evolution of function dx<f for the same set of parameters as in Figure 
but for an increasing L{t). The initial perturbation is damped (Left), but some 
instability arises as L{t) gets larger after some time (Right). 
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Figure 7: Evolution of function dx^p for an unstable set of parameters a = 5.10“^ 
and d = 4, and two different initial conditions. Final shape plotted in dash line 
on the left corresponds to the surface plotted on the right. The color (online) 
represents the quantity of the density of materials available on the membrane. 


observe clearly the cigar shape up in Figure and the rod shape below in 
Figure The latter is typical of fission yeast growth. 
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A Kinematics 


In this appendix, we give more details on the way to derive equation (2.2) 


governing the evolution of the generatrix curve Ct pushed by a vector fields, 
v(t, •) : Cf 


Let t be a positive real and s E [0, L{t)], we express ^ 


d / dr{t,s{t))\ _ 
dt\ ds ) ~ 

1 /r(t + dt, s + d§) — r(t + dt, 5) dr{t,s)\ 
dt,ds^o dt V ds ds ) ' 


(A.l) 


where dt (resp. ds and then ds) is an infinitesimal time (resp. length). 
We denote by s the curvilinear abscissa at time t and by s the one at time 
t + dt. 

However, we have (see Figure]^ 


r{t + dt, s) = r{t, s) + (v • ei) dt + o(dt) , 
r{t + dt, s + ds) = r{t, s + ds) + ((v + dv) • ei) dt + o(dt), (A.2) 

here the vectors dtv (resp. dt (v + dv)) represents the vector fields v(t, •) 
at (r(t, s), z{t, s))) (resp. at (r(t, s + ds), z{t, s + ds))) as depicted in Figure 

El 


Using the relations (A.2), we get 


r(t + dt, s + ds) — r(t + dt, s) fdr{t,s) (dv-ei)dt\ds 


ds 


- 


ds 


+ 


ds 


ds 


+ o(dt). 


(A.3) 


Notice that the derivative of the vector field v with respect to 

d'v{t,s) (dv-ei) 

ds ds 

Moreover, the infinitesimal quantity ds reads (see Fig ) 


s is 


(A.4) 


ds = ds — dt (v • r) + dt ((v + dv) • r) + o(dt) , 
= ds + dtdv • T + o(dt) 


By dividing the previous equality by ds and inverting it, we get 


— - 1 

ds ds 


T + o(dt). 


(A.5) 
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By plugging (A.4) and (A.5) into (A.3[), we obtain 


r(t + dt, s + ds) — r(t + dt, s) dr{t,s) d'v{t,s) 

+ dt-- • Gi 


ds 


— dt 


ds ds 

dr{t, s) d'v{t, s) 


ds 


ds 


T + o(dt) . (A.6) 


By plugging (A.61 into (A.l | and by expanding ^ the 

equation, governing the evolution of the generatrix curve for all t > 0 and 
for all s G [0, L{t)], 


d fdr{t,s)\ dsd‘^r{t,s) dr{t,s) fd'v{t,s) 


dt 


V ds ) dt d'^s ds V 

.dv{t,s’) 


dt d'^s 
ds 


ds 


r(t,s)^ + 


(9v(t, s) 
ds 


•ei . 


ds 5 ) f f 

Finally, noticing — = / ( - -4—) • r(t, s') ds', we obtain for all t > 0 

dt JO ds _ 


and for all s G [0,L(t)] the equation (2.2). 


B Change of variables 

This section is devoted to some preliminary simplifications of the system 


of equations (2.2|-(2.6)-(2.8|. Firstly, we change the parametrization of the 


generatrix curve and derive its associated evolution equation. This leads to 


simplify the evolution equation (2.2). Furthermore, through changing vari¬ 


able, we bring back the study of the system (2.2|-(2.6|-(2.8) defined on a 


varying-time domain into a fixed one. Finally, we combine the velocity vec¬ 


tor fields equations (2.6) with the evolution one to get a condensed formula 


for the system of equations (2.2 )-(2.6)-(2.8) 


Change of parametrization 

In what follows, instead to reconstruct the generatrix curve by using the 
system of coordinates {(s, r(t, s))}, we will use the parametrization given by 
{(s,V9(t,s))}. 


By using relation (3.1), equation (2.2) can be expressed as a partial dif¬ 


ferential equation on (t, s) i—)• ip{t, s). Indeed, the left hand side of equation 


(2.2) reads, 

Vt > 0 , Vs G [0, L(t)] , dt {dsr{t, s)) = — sin((/;(t, s)) dt(p{t, s). (B.l) 

By using the expression of the two vector fields n and r (given in Fig 
[^, we get 


V • ei = sin((^) -|- Vr cos{ip) , \r ■ r = Vn dsip + dsv-r 


(B.2) 
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Consequently, by plugging the relations ( |B.2 ) into the right hand side of 
equation (2.2) we obtain Vt > 0,Vs G [0,L(t)] , 


ds yj (9sv) • Tj dsr + v • ei j = 

sin((/7) ^ J Vndsi^ dsip — dsvj^ 


(B.3) 


Then, by dividing (B.l) and (B.3) by sin(</5), we have for all t > 0 and 
s G [0, L{t)] 

dt<p{t, s) = (^- Vn{t, s)ds(f{t, s')ds'^ ds(p{t, s) - dsvnit, s) , (B.4) 

where Vn satisfies the equation (2.6bI. 


Moreover, it derives from relation (3.1) that for all t > 0 and for all 
s G [0, L{t)], 


r{t,s)= / cos((/9(t, s'))dsC 

Jo 


Ks{t, s) = dsif , 


Keit,s) = 


sin((/?) 


cos{(p{t, s'))ds' 


By using relations (B.5), equation (2.7) becomes 


(B.5) 


/3 


/ \ 

dsip{t, s) sm{ip{t, s)) 


V 


cos{(p{t, s'))ds' 


5s cos{ip{t,s'))dsdsfj,{t,s] 


cos{ip{t, s'))ds' 


+ s) 


(B.6) 


Change of variables 

The system governing the dynamics of the cell membrane is posed on a 

S i^—I i^—I 

time-varying domain [0,L(t)]. By setting x := the system (B.4)-(B.6) 


can be expressed on the fixed domain [0,1]. 


L{ty 


In what follows, for all the function / defined on M'*' x [0, L(t)], we denote 
by / the function defines on M+ x [0,1] such that 

Vt>0, VxG[0, 1], f{t,x) = f{t,xL{t)). 
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By replacing s by xL{t) and using Definition ??, the equation (B.4| 
becomes for all t > 0 and Vx G [0,1], 


n ^ fdMt) (r~u ~u 

dMt,x) =1 \Jq '^n{t,x)dxip{t,x)dx j — 

dxVn{t,x) 

L{t) ’ 

As well, the normal velocity field Vn is equal to 

Vn{t, x) = COt{ip{t, x))vt{t, x) + H{t, x) , 
where the function H is defined on M"*" x [0,1] as 

H{t,x) = {de{t,x) - uds{t,x)) . 

K0{t,x) 


(B.7) 


(B.8) 


Moreover, since Vr is solution to the ODE (2.6b), the tangential velocity 
field vt is the solution of the following ODE 

dxVt{t, x) - (^dx x)) cot((f(t, x)) j vt{t, x) = d{t, x) (B.9) 
where the function d is defined on M'*' x [0,1] as 

d{t, x) = L{t)^{t, x) [ ds{t, x) — vde{t, x) ) — L{t)Rs{t, x)H{t, x). 


Let us expand, by using the equalities (B.8) and (B.9), 


px 

/ Vn{t,x)dx(p{t,x')dx = 

Jo 


-Vt + \^d{t,x) + H{t,x)dx{(p{t,x))jdx 

However, by differentiating Vn with respect to x, we get 
dxVn = dx(p Vt - cot{(f) d + dxH 


(B.IO) 


(B.ll) 


Then, plugging (B.lO)-(B.ll) into (B.7), the function (p is solution to 


the following partial differential equation 
'dtL{t) 




+ cot{(p{t, x)) d{t, s) — dxH{t, x). 


(B.12) 


Let us call J1 = Finally, equation (B.6) could be expressed for all 

times t on a fixed domain [0,1] as 
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with the boundary conditions lcos{t,0)dxfl{t,0) = Icos(t, 1) = 0 

2 

and where a{t) = . 
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